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the  product  of  a distance-graph  with  a flow-graph  several  special  cases  of  [ 
the  QAP  are  Investigated.  A polynomial -growth  algorithm  is  desv.ribed  for  [ 

the  QAP  when  the  distance  and  flow-graphs  are  Isomorphic  trees.  In  tlie  case  i 
when  the  graphs  are  single  stars  the  algorithm  becomes  the  well  known  rule 
for  multiplying  two  sequences  of  numbers.  The  case  of  a complete  distance- 
graph  and  a tree  flow-graph  becomes  the  travelling  salesman  problem  when  the 
tree  is  a hamiltonian  chain  and  the  flows  are  all  unity.  A dynamic  programming 
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of  mi  nimi  7.  the  product  of  a distance-graph  with  a flow-graph  several 

special  cases  of  tlie  QAP  are  investigated.  A polynomial -growth  algorithm 
is  described  for  the  QAP  when  the  distance  and  flow-graphs  are  isomorphic 
irees.  In  the  case  when  the  graphs  are  single  stars  the  algorithm  becomes 
the  wi.'ll  known  rule  for  multiplying  two  se^iuences  of  numbers.  The  case 
III  a coni|ilelc  distance-graph  and  r.  tree  flow-graph  becomes  the  Iravellin:' 
salesni.in  prohleiii  when  the  tree  is  n hamiltoninn  chain  and  the  I lows  are 
,1 ! I nmtv.  A dynamic  programming  algorithm  is  presented  for  tlu'  case  oi 
the  I low-graph  being  a general  tree  with  arbitrary  flows.  The  very  special 
case  of  "narrow"  bipartite  graphs  is  also  considered. 


1 . ] NTRODUCTION 

Consider  p machines  1 , . . . .a, . . . , p with  a known  flow  of  material 

f ^ between  every  pair  of  machines  (a, 3).  Let  there  be  q ■■  p locations 
OP  ~ 

1 , . . . , i , . . . , q ith  known  distances  d^^  between  every  pair  of  locations 
(i,j).  An  assignment  of  machines  to  locations  is  a one  to  one  mapping 
0 of  the  set  of  machines  into  the  set  of  locations,  so  that  p(g-)  is  the 
location  that  machine  a is  assigned  to. 

Tlic  cosl  of  a mapfilrig  p is  <leflned  as 


zip')  ■-  ^ f 


0-3  pialpih) 


(liven  the  two  matrices  [f  | and  fd..l,  the  quadratic  assignment 

Q'P  1 j 

problem  (QAPl  is  that  of  finding  a mapping  P*  which  minimizes  z(p)  as  given 


by  (11 


The  QAP  appears  in  a number  of  special  location  problems  such  as 


I lie  allocation  of  machines  to  locations  - used  above  to  introduce  the  QAP  - 
the  location  ol  electronic  components  on  cinnit  boards  [1?1,  the  ordering 
oi  I nt  errela  teil  data  on  magnetic  tap<»,  etc.  Other  examples  not  invnlvin}'. 
special  location,  but  which  can  be  lormulatcd  as  QAP's  include  the  trian- 
gulai  ization  of  economic  i nput -outpu  t matrices  fl’l,  the  minimization  of 
average  job  completion  time  in  machine  scheduling  [8]  and  extensions  of 
the  travelling  salesman  problem  [7], 

A survey  of  exact  algorithms  for  the  general  QAP  is  given  by  Pierce 
and  ( rowston  [11],  and  an  improved  algorithm  i..  described  ir  f‘^1-  Exact 


algorithms,  however,  are  unable  to  solve  general  QAP's  of  even  moderate 
size  [4].  Approximate  algorithms  for  the  QAP  are  surveyed  by  Nugent  et  al 


[10|  and  Moore  [4],  while  Sciabin  and  Vergin  [13]  demonstrate  that  these 
are,  in  general,  unsatisfactory. 


JT  J r-7 ! ;^  f:// "7/j/- 


In  till';  p«|i(M  w(‘  lonsliler  special  cases  of  the  (^AP  which  /ire  e/isler 


In  solve.  lU'  iecaslin>;  the  ()AP  in  graph  theoretic  terms  as  Llie  dot-prodiict 


>1  a di SI ancc-graph  with  a flow-graph,  we  consider  cases  when  these  graphs 


have  special  I onus . In  particular  we  describe  a polynomial  growth  algorithm 


lot  the  QAP  when  both  of  these  graphs  are  trees.  Vfhen  only  one  of  the  two 


graphs  is  a tree  and  the  other  is  a complete  graph,  the  QAP  can  be  solved 


!-v  a special  dynamic  programming  algorithm  which  is  a generalization  of  a 


similar  algorithm  for  the  travelling  salesman  problem.  This  last  case  occurs 


veiv  often  in  [iractical  location  prob lems , e . g . , in  the  layout  of  an  assembly 


-3- 


2.  GRAPHICAL  REPRESENTATION 


A graph  G is  defined  by  the  doublet  (X,A)  where  X is  a set  of 


vertices  and  A a set  of  links.  Unless  otherwise  specified  we  will  use 


’graph"  to  mean  a "non-directed  graph  without  loops."  The  terminology 


used  is  I roin  [ 1 , 1]  . 


Given  a grapli  G*  (x'.A^l  with  a cost  matrix  isomorphic- 


graph  C.  " - (X'  , a 'I  with  a cost  matrix  ® mapping  p of  X onto 


the  dot -product  graph  is  written  - using  the  product  operator  rifpl  - as; 


G nCplG 


and  is  defined  as  the  graph  G = (X,A)  isomorphic  to  G with  costs  given 


bye  . . . 

1!  i]  p(i)p(.ll 


The  va lue  of  a graph  G - (X,A>  is  defined  as: 


V(G1  - 


( X , , X . ) eA 
i ,I 


An  image  of  a graph  (i  ^ in  a graph  (i  is  anv  partial  subgraph  of  G^  which 


is  isomorphic  to  g\  Wg  will  denote  by  M(G  \g  the  family  of  all  sucli 
image  graphs.  The  cardinality  of  the  set  M(G  \g  is  called  the  image 


number  of  G^  in  G and  is  denoted  by  m(G  ,G  ^ ). 


The  QAP  can  now  be  restated  in  the  following  way. 


Let  G^  = (X^.A^l  be  a flow-graph,  whose  vertices  X^  represent  the 


set  of  machines  and  the  link  costs  are  the  flows  between  the  corresponding 


machines.  Similarly,  let  G^  = (X*^,a'^')  be  a di.--  ance-graph,  whose  vertices 


represent  tlie  set  of  locations  and  the  link  costs  are  the  dis  ance;.  between 


the  cr rresponding  locations.  Wo  will  assume  (without  loss  of  generality) 


that  jx**!  |X^ 


n 


H 


H 
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The  QAP  is  Chen  the  problem  of  finding  a graph  G and  a mapping 


p ol  G on  G which  minimizes  the  expression; 


Mi  II  I Min  V(C  n(pU:)  i 


I .1  P 
(G:M(G  ,C  ■> 


Tlte  niiinlitT  nl  il  I f I ereni  mappings  p nf  G onto  an  Isomorphic  graph  is  the 


isomorphic  mimher  sic^)  of  G^.  Thus,  the  inner  minimization  of  (,'M  is 


nvor  n set  of  cardinality  s(G  ) and  the  outer  minimization  is  over  a set 


f d 

of  cardinality  m(G  ,G  ) 


It  is  reasonable  to  expect  that  as  s(G^')  and  m(G^,G*^')  increase, 


the  difficulty  of  the  QAP  will  also  increase.  Cases  when  both  s(G^'l  and 


ml(;  ,G  ) are  small  can  be  solved  trivially  bv  enumeration.  Very  few 


i|H'cial  cases  in  v-?hich  onlv  one  of  these  two  numbers  is  large  can  be 


solv'd!  Iiv  polvnomlal  growtii  a I g.or  i i hms  . The  case  whert'  (/  (i'*  K f the 


vomplele  grapli  on  n vertices)  is  the  problem  usually  considereil  in  tin 


iterature  as  the  general  QAP  ami  has  miG^G.'*)  - ] and  s(G^)  = n! 


i.  CASKS  WITH  C K . (Imago  niimbor  - 1) 


WluMi  I lio  Imago  minil'or  in((l  ,C  ) 1 tlio  outer  minimization  ol  ( ,’ ) 


beeomos  roilundant'  and  only  the  inner  minimization  remains. 


1.!.  Trivial  cases  (small  isomorphic  number) 


Trivial  cases  that  can  be  solved  by  complete  enumeration  are; 


(i)  Chains : When  C.^  and  are  chains,  the  isomorphic  number 


s(G  i = 2 and  the  inner  minimization  in  expression  (?)  onlv 


involves  two  evaluations. 


(ii)  (!yc  les : The  isomorphic  number  s((;  ) = ?n. 


(iiil  Whce  I s : The  isomorphic  nuinoer  s(G  ) = 2(n-l). 


(iv)  Regular  graphs:  Certain  regular  graphs  (e.g.,  webs  of  low 


order'/  have  small  isomorphic  numbers  and  can  be  enumerated. 


There  is,  however,  great  variation  in  the  isomorphii.  numbers 


of  regular  graphs  even  of  the  same  degree  as  shown  bv  the 


example  in  Fig.  1.,  and  no  general  statement  can  be  made. 


i.2.  Solvable  cases  (large  isomorphic  number 


A.  Simple  stars:  Wlicn  C*  and  (/*  are  simple  stars  with  one  central 


vertex  (with  index  0)  and  n outer  vertices,  the  isomorphic  number 


is  n!  However,  this  special  QAP  can  be  solved  by  a well  known 


rule  namely:  Order  the  n flows  f in  ascending  order,  and  the 

Oq 


n distances  d^^  in  ilescending  order.  The  optimum  mapping  p* 


then  maps  the  flow  in  the  flow  list  to  the  distance 


in  the  distance  list  for  all  k = l....,n. 


B.  Multiple  stars 


The  graph  in  Fig.  2 shows  a 28-vertex  3rd  order  star 


with  vertex  1 as  the  center.  Consider  a general  k-order  multiple 


...■'I 


W’  -■ 


TvaJO  R 


h(  I'.oMMk  [->1  \C 


CjRaphs  or 


ami  let  V . be  the  value  of  the  solution  of  this  problem 


Update  c(a,i>  = f , ^ . d ....  + V . 

p(o)o  P(i)i  Ql 


tiYWrWiiriiil 
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star  rrapli  arbitrarily  rooted  at  the.  center  vertex  and  let  the 
l;il)(*l  ( X I III  vertex  x lie  the  cardinality  of  th(*  path  from 
(III-  iiiilei  Id  V.  Mil-  l/ib(‘l  of  tile  Dul(*rmo.ML  vtirtlcf.'i  In  then 
k and  the  label  o!  the  center  vertex  is  zero. 

We  describe  below  a dynamic  programming  algorithm  for  the 


solution  of  QAP's  involving  arbitrary  k-order  stars. 


Let  the  flow  and  distance  graphs  be  = (X^,A^)  and 


t;^  fx'^.A^’)  respectively.  For  any  vertex  x e and  x^  e X*^ 


let  clj,!")  be  the  minimum  cost  of  mapping  x and  all  its  successors 

Q 


(i.e.,  vertices  reachable  from  x via  arcs  of  the  rooted  tree) 

Q 


to  X.  and  all  its  successors.  We  will  denote  by  x , , the 

1 p(a) 


predecessor  ol  vertex  x 


Description  of  the  a I >',or  i thm  (lor  k-oriler  stars) 


Step  1.  For  each  >;  e X*  and  x.  e x'^  with  -(.(x  ) = (x  ) = k set 

0 i O'  i 


r,  , . d , . , . 

p(qV<  p(i)i 


Set  LEVEL  - k - 1 


Step  2.  For  each  x c and  x.  e X^  with  ) = -tCx.)  = LEVEL  calculate 

a 1 O’  1 


c(o,i)  as  follows: 


Let  --  iPjpO)  = a] 

and  = '.ilp'T'  ^ i] 

lii)  Set  up  the  linear  assignment  problem  with  cost  matrix 


M(i  . r hikii  i.i  . \ A w..  viw  w tvl " 

1 C.0NUIK)"\WC  k Z^:.6M  4SS 


/(  lnt  t:R 


■<X 


}(  ' o 


/ 


x>" 
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Step  3.  Set  LEVEL  = LEVEL  -1.  If  LEVEL  = 0 go  to  (4)  else  go  to  (2). 

Step  4.  Stop.  If  X and  x.  are  the  center  vertices  of  and 

o 1 

o o 

respectively,  Is  the  value  of  the  solution  to  the  QAP. 

(Note : Tlie  mapping  corresponding  to  this  solution  value  can  be 

found  hv  backtracking  in  the  usual  dynamic  programming  manner.) 

I'he  above  algoritiim  is  good,  exhibiting  polynomial  rate  of  growth 

with  tiie  total  number  oi  vertices  in  the  k-ortier  star.  Thus,  if  each 

k+1 

vertex  of  the  star  has  exactly  m successors,  there  are  (m  -l)/(m-l) 

k-1  2k-2 

vertices  in  all.  The  algorithm  involves  2m  sortings,  m evaluations 

2k-'>  2 

and  the  solution  of  (m  -l)/(m  -1)  assignment  problems  of  si7.e  m X m. 
ft  shouUi  be  noted  here  that  in  the  first  pass  through  step  2 of  the  above 
■ilgr)rithm,  the  solution  of  the  assignment  problems  defined  by  C is  unnecessary 
since  these  assignment  pi-oblems  correspond  to  simple  stars  and  do 

solved  by  the  simple  ordering  described  earlier.  Thus,  in  the  case  of  the 
simple  star  the  above  algoritiim  disintegrates  to  the  well  known  ordering 
ru  le  of  3 , ’ ( A ) . 

C . Ceneril  trees 

Any  tree  can  be  arbitrarilv  rooted  and  considered  as  a k-order 

star  with  «>-cost  ( f low/d  ist  ance)  arcs.  The  algorithm  described  above 

can  therefore  be  used  to  solve  i).-\l’'s  with  general  trees.  This  ir  equivalent 

to  slightly  mollifying  the  above  algorithm  so  that  c(q,1)  is  finite  only 

for  those  pairs  of  vertices  [x  ,x.|  for  which  the  subtree  defined  by  x 

o 1 a 

and  its  successors  in  G is  isomorphic  to  the  . ibtree  defined  by  x^  and 

. 

Its  successors  in  G . 
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n.  Narrow  Bipartite  Graphs 

Cons1(l(>r  tlic  (^Al’  when  ami  (1^  are  complete  bipartite  graphs, 
say  of  the  form  K . Using  the  normal  notation  we  will  express  in  terms 

1 j 

of  its  two  independent  vertex  sets,  i,e.,  = (X^,X^)  and  similarly 

= (x'^.x'^).  We  call  K narrow  if  min  (r,s)  ■'  max  (r.s).  Let  us 

r s r , s 

assume  that  r - min  (r,s)  and  specifically  consider  the  case  when  r is  small. 

The  isomorphic  nuinher  of  K is  r!sl,  however,  if  r is  small  enough,  the  r! 

r , s 

possible  mappings  of  X^  on  X^  could  be  enumerated.  For  each  such  mapping  o 

we  would  then  compute  ¥ x e X^  and  x,  e X*^: 

a « IS 

ii  r 

and  solve  the  s by  s linear  assignment  problem  with  cost  matrix  [c(a,i)]. 

Tlie  least  cost  assignment  solution  over  all  r!  mappings  p i then  ll.'  solution 
to  the  QAP. 


Obviously,  such  a procedure  is  only  j>racticnl  when  r is  very  small 

V 5 

(say  _ ll  l)ut  with  a g.iven  r t lie  complexity  as  a Tunction  of  s is  Ois'  ' ) 


since  it  only  involves  the  linear  assignment  problem. 


''(.  CASI’S  WIIHN  C IS  COMPLKTK 


W('  w(  I 1 now  I .'lUc  C Ic)  lit'  ;i  cninplrti-  ^rMph  on  ii  vi‘rCic(‘H  /in<l 


Ciinsldor  cases  vt/hen  C.  Is  ol  tllllt’rc'nt  fomia. 


^'  _ iJLlL  simple  star  with  n vertices 


The  iinane  number  m({;  ,(l  ) is  n.  I'^ach  Imaj'e  corresponds  to  a 


star  partial  graph  of  G with  a specific  center  vertex.  Once  an  image  graph 


is  cltosen  the  optimum  mapping  p of  onto  G^  can  be  found  as  in  section 


3.2(A)  earlier.  The  total  complexity  of  the  procedure  is  therefore  0(n  log  n) 


B.  General  trees  with  n vertices 


Altliough  the  procedure  for  simple  stars  given  in  section  3 for  the 


case  in  whicli  G*^  - G^  is  generalized  above  to  the  case  in  which  is  a 


complete  graph,  the  corresponding  algorithm  of  section  3.2(B)  for  k-order 


stars  (or  arbitrary  trees)  does  not  g.enernlize.  llie  fact  that  sucli  .1  generaliza- 


tion is  not  possible  can  1)C  demonstrated  by  considering  G to  be  a simple 


chain  of  n vertices  numbered  consi'cuti vely  I rom  an  end  vertex  and  take  f 


0 ,crH 


for  all  a ~ l,...,n-l.  We  now  have  s(G*)  - 2 and  m(G*^,G*^)  = l/2n!  In  fact, 


the  image  graphs  of  G^  in  G*^  are  all  the  hamiltonian  paths  of  g'^,  and  since 


we  have  taken  all  flows  to  be  unity,  the  value  of  the  product  graph  of  G 


with  an  image  graph  is  simply  the  length  of  the  hamiltonian  path  forming  the 


image  graph.  (Note  that  the  2 possible  mappings  of  G onto  the  image  graph 


give  the  same  value.)  Thus,  the  QAP  with  G a simple  chain  and  G the  complete 


graph  becomes  ofiuivalent  to  the  open-ended  travelling  salesman  problem. 


Althoug’h  tbe  algorithm  of  section  3.:’(B)  do  s not  generalize  to  the 


present  case,  this  case  (of  ('.*  being  an  arbitrary  tree  and  G^*  tlie  comiilete 


graph)  is  possibly  the  most  important  of  all  cases  of  the  QAP  - as  far  as 


practical  applications  are  concerned  - since  in  many  situations  (e.g,  , in 
assembly  line  layout,  pipeline  design,  etc,)  the  flow  graph  is  of  this 
form.  In  view  of  its  importance  we  present  here  a specialized  algorithm 
which  can  solve  QAP' s of  considerably  larger  size  than  any  algorithm  for  the 
general  QAP. 

The  algorithm  is  a generalization  of  the  dynamic  programming  algorithm 

in  [6]  for  the  travelling  salesman  problem.  The  generalization  is  in  two 

directions:  (i)  it  considers  different  flows  between  machines,  and  (ii)  it 

can  accommodatf  arliitrary  trees  instead  of  simply  chains. 

Consider  a general  tree  graph  and  suppose  it  is  arbitrarily 

rooted  at  some  vertex  x , Let  T(x  ) be  the  directed  subtree  reachable  from 

o a 

X , including  x itself  as  the  root  of  the  subtree.  We  will  also  use  T(x  ) 
a a a 

to  mean  the  set'  of  vertices  of  this  subtree.  Let  x , , be  the  immediate 

P(a/ 

predecessor  of  vertex  x in  the  rooted  tree  G . 

O' 

For  a given  vertex  x e let  S C X*^  and  x.  e S where  Is  I = 1t(x  ) 1 . 

Q a X O'  ' o'  ' O'  ' 

Let  ■ S ’be  rhe  induced  subgraph  on  the  subset  S of  vertices  of  the  distance 
O'  O' 

graph  . Define  the  function  g(S  ,x.)  as: 


R(S^,x.) 


Min 

C€M(T(x  ),  S ■•) 


Min  V (T(x  )n(p)G) 

P 

-p(x  )=x. 

O'  1 


i.e.,  g(S  ,x.)  is  the  solution  to  the  QAP  defined  on  the  subgraphs  T(x  ) and 

O'  X O' 

• ' S with  the  restriction  that  the  optimum  mapping  should  have  p(x  ) = x.  . 

O'  O'  X 

The  function  g(S  ,x.)  can  be  computed  recursively  as  follows: 
a X 

(i)  If  X is  the  predecessor  of  only  one  vertex  Xg.- 

a 

g(S  ,x  ) = Min  [g(S  ,x  ) + f d 1 (4: 

ox  P j OP  ij 

j P 

where  S = S.  U [x.  } . 
a p ^ x-* 


1 

if 

1 -1?- 

1 

1 

^ f 

i.  QAP  with  vertex  1 of  the  flow  f>raph  G mapped  onto  vertex  x of  the  distance 

1 

1 d 

f graph  G . The  value  of  the  solution  to  the  QAP  would  then  be: 

1 

1 

FI 

’ z = Min  [g(S  ,x  )] 

r „d  1 1 

1?  x.cX 

■ 1 

'1 

1 f 

1 It  is  interestir  to  note  that  the  case  of  G being  a single  star- 

F graph  involves  only  one  application  of  iteration  (5)  to  calculate  g(S^,x.) 

1 

!:  for  eacli  x and  oven  for  that  single  application,  the  linear  assignment  problem 

e ^ 

>3 

can  be  used  to  solve  the  outside  minimization  as  mentioned  earlier.  This  is  the 

■J 

1 

1 simplest  case  mentioned  in  section  4. A.  above.  The  computationally  most  difficult 

■i  f 

1 case  is  the  case  of  G being  a hamiltonian  chain  which,  as  mentioned  earlier, 

K 

; leads  to  the  travelling  salesman  problem.  Cases  of  trees  G with  values  of 

1 graph  diameter  1 etween  these  two  extremes  are  of  intermediate  complexity. 

1 

r C,  Narrow  bipartite  graphs 

F f 

1 When  G = (X  ,X  ) is  the  complete  bipartite  graph  K with  r ■ s the 

X.  r s IT  ^ s 

? 

f. 

i QAP  can  he  solve  d when  r is  very  small  by  the  method  of  section  3.D.  , i.e.. 

1 

b simply  enumerating  all  r!(^)  mappings  of  the  set  X*^  into  a set  of  r vertices 

‘J-  Q 

F of  G and  solving  an  s by  s assignment  problem  for  each  such  mapping.  Obviously 

1 

1 this  would  only  be  practical  for  r = 2 or  at  most  3 even  for  graphs  G^  with 

1 

1 only  20  or  30  vertices, 

f f 

D.  (.'  is  a collection  ot  links 

1 

: j- 

1 For  Instructional  purposes  it  may  he  worthwhile  to  note  that  if  G 

:i 

K 

f is  a disconnected  graph  composed  ol  q components  each  of  which  is  a single 

1 link,  then  if  all  flows  are  unity  the  QAP  becomes  a matching  pri.hlem  and 

if 

1'  can  be  solved  as  such.  It  is  not  at  all  clear  if  the  problem  with  non-unitv 

ji 

i 

^ 1 

1 ■ 

If 

S' 

li 

1 

n 

^11  1 . , 

w m 
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flows  could  be  solved  as  a matching  problem.  Although  the  matching  algorithm 


can  also  be  used  to  solve  general  degree-constrained  partial  graph  problem,^ 


defined  on  c‘' , the  form  of  the  sol.it  ion  cannot  be  guaranteed  to  correspon.l 


to  any  a priori  defined  flow  graph  G . 


V ... 


5.  CONCLUSION 


We  have  ex(n'essecl  the  f)AP  in  tenas  of  Rraph  multiplication  and 


clasfUfied  and  invest  i p.at  ed  c.ises  d(»pendinR  on  the  form  that  the  distance 


graph  c'*  and  Flow  graph  C*  Lake.  Tlie  case  when  C^  and  C*'^  are  both  simple 


st.ir  graplis  was  known  as  ;i  solval»l(‘  case  of  the  QAP.  A pol ynomia  1 -growth 


algorithm  has  now  been  given  for  the  solution  of  QAP's  when  both  C^  and  C^' 


are  arbitrary  tree.s.  Although  the  algorithm  for  simple  star  graphs  generalizes 


to  the  case  of  one  graph  being  a star  and  the  other  being  a complete 


graph,  the  new  algorithm  for  arbitrary  trees  does  not,  since  any  such  generaliza- 


tion implies  the  travelling  salesman  problem.  However,  a specialized  algorithm 


for  solving  QAP's  where  G is  any  tree  and  G the  complete  graph,  is  described 


which  can  solve  considerably  larger  problems  than  any  general  QAP  algorithm. 
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